Emergence of turbulence in an oscillating Bose-Einstein condensate 
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We report on the experimental observation of vortex tangles in an atomic BEC of ^'^Rb atoms 
when an external oscillatory perturbation is introduced in the trap. The vortex tangle configuration 
is a signature of the presence of a turbulent regime in the cloud. We also show that this turbu- 
lent cloud has suppression of the aspect ratio inversion typically observed in quantum degenerate 
bosonic gases during free expansion. Instead, the cloud expands keeping the ratio between their 
axis constant. Turbulence in atomic superfluids may constitute an alternative system to investigate 
decay mechanisms as well as to test fundamental theoretical aspects in this field. 
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Turbulence in classical fluids has always been an in- 
triguing issue in physics and understanding and control- 
ling turbulence represents one of the main challenges in 
physics p. Nevertheless, classical turbulence is a con- 
cept difficult to define and there are no fully consistent 
theories to describe it A turbulent flow is a highly 
unsteady flow characterized by the presence of vorticity, 
diffusion and energy dissipation processes. In a turbulent 
fluid, the fluctuation of pressure and velocity occur in a 
very broad range of length and time scales, making its 
mathematical description very complex 

Turbulence is also possible in superfluids, such as su- 
perfluid Helium 4] and superfluid ■^He-B Turbulence 
in superfluids or Quantum Turbulence (QT) is character- 
ized by the appearance of quantized vortices distributed 
in a tangled way 0, 0, not forming regular lattices. 
This subject has been extensively studied in superfluid 
Helium both experimentally and theoretically since it has 
been discovered about 50 years ago Until recently, 
turbulence in He-superfluid could only be observed by 
indirect methods. Oscillating wires or wire grids [lO| 
were the most commonly employed techniques. Recently, 
there was the first direct observation of QT in Helium 
llj . That work opened upnew frontiers on the study of 
QT in superfluid Helium [l^]. In the theoretical branch, 
QT in Helium has been satisfactorily modeled by the 
vortex filament model [l^] , though some features such as 
vortex reconncctions must be introduced artificially. 



Nevertheless, the achievement of BEC in trapped 
atomic gases [l4[ and the subse quen t observation of quan- 
tized vortices in these samples [l^ opened up the possi- 
bility to study turbulence in a more controlled fashion, 
since one can control the main characteristics (interac- 
tion, atomic density, number of atoms, trapping config- 
uration) of the atomic sample. Additionally, the possi- 
bility of expanding the quantum cloud by releasing the 
atoms from the trap, makes possible the direct observa- 
tion of the vortices by optical means. There are today 



many open questions related to QT, including mecha- 
nisms of decay Q which may find important hints on the 
direct observation of QT in a trapped atomic superfluid. 
In fact, the study of QT in atomic superfluids may shine 
light on the turbulence characteristics that are universal, 
due to the superfluidity as a macroscopic quantum state, 
and reveal which characteristics depend on the specific 
nature of the superfluid. Specifically, in atomic superflu- 
ids, the ratio of the average vortex spacing to the vortex 
core size is expected to be much smaller than in Helium, 
leading to different modes of turbulence decay. Several 
extensive reviews 13, 16, 17| summarize the present sta- 
tus and challenges to be faced in the research field of 
QT. 

Shortly after trapped BECs were first realized, the ob- 
servation of vortex lattice structures, as well as the crys- 
tallization dynamics of these structures were reported 
(ist - These dynamics have been successfully confirmed 
quantitatively by the numerical simulation of the Gross- 
Pitaevskii equation [19|. However, in the experimental 
research of trapped BECs, turbulence has not yet been 
investigated. Recently, this topic has attracted theoret- 
ical interest 0, IIqI, HH and numerical simulations show 
that turbulence is probably present as a step in the for- 
mation dynamics of a BEC with the presence of spon- 
taneously formed vortices 2M- Nevertheless, no experi- 



mental evidence of this regime has been specifically pre- 
sented so far. In this Letter, we report on the experimen- 
tal observation of QT in a magnetically trapped BEC of 
^^Rb atoms evidenced by the observation of tangled vor- 
tices in the quantum sample, an intrinsic characteristic 
of QT, as well as a change in the hydrodynamic behavior 
of the atomic cloud. 

The experimental sequence to produce the BEC as well 
as the procedure to generate vortices in the condensate 
are described in details in Refs.[23l. [23 |. In brief, we 
produce a BEC of ^''Rb containing (1 — 2) x 10^ atoms 
with a small thermal fraction in a cigar-shaped magnetic 
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trap with frequencies given by = 27r x 210Hz and 
LOx = 2tt X 23Hz. Once the condensate is obtained and 
still held in the trap, an oscillating magnetic field is ap- 
plied. This field is produced by a pair of anti-Hclmholtz 
coils placed with their axis closely, though not exactly, 
parallel to the weak axis of the trap (FiglT^). Since the 
coils are not perfectly aligned to the condensate axis, the 
external field in the region of the trap has components 
parallel to each trap eingenaxis. Additionally, the com- 
ponents along the two equal directions that give the ra- 
dial symmetry of the trap are slightly different, breaking 
it up as well as the symmetry of oscillation of the conden- 
sate. This oscillatory scheme was developed to investi- 
gate coherent mode excitations in Bose-condensates [i^ . 
Such excitations are also believed to produce vortices. 

A pictorial view of the field acting on the condensate is 
shown in FigUlD. A combination of shape modification, 
displacement of the minimum and rotation is coupled to 
the atoms. The quadrupolar field is applied for a period 
ranging from 20 to 60 ms, oscillating sinusoidally in time 
with a frequency of 200Hz. This extra magnetic field has 
an offset so that in a period of oscillation, it goes from 
zero to a maximum value and back to zero, never invert- 
ing its direction. The maximum value reached by the 
field gradient characterizes the strength of the excitation 
and has been varied from zero to 190 mG/cm in the axis 
of the anti-Helmholtz coils. After the end of the oscilla- 
tion stage, the atoms are left trapped for an extra 20 ms 
before being released and observed in free expansion by 
a standard absorption imaging technique. For small am- 
plitudes of the oscillating field as well as short excitation 
periods we observe dipolar modes, quadrupolar modes 
and scissors modes 2^ 27 1 of the BEG but no appear- 
ance of vortices. Increasing both parameters we start to 
see vortices that grow in number with amplitude and/or 
duration of excitation. These regimes are well described 
in a previous report [il ]. 

The effective mechanism that nucleates the vortices is 
still under theoretical investigation [2^. We believe that 
a possible mechanism of vortex formation are the so- 
called Kelvin-Helmholtz instabilities [2^, a well-known 
phenomena in the scope of fluid interface theory and ex- 
periments. These instabilities occur at the interface of 
two fluids with a relative motion and give rise to the for- 
mation of vortices. This phenomenon has been theoreti- 
cally investigated [30| in normal fluid-superfluid mixtures 
and directly observed in mixtures of A and B phases of 
superfluid liquid Helium [3l|. The fact that this mech- 
anism may be taking place is supported by previo usly 
theoretical and experimental observations 21,24, 2^, 27 1 . 
Ref.[2l| shows that a stirred condensate ejects energetic 
atoms to form an outer cloud durin g the q uadrupolar 
mode excitation. Additionally, Refs. 2^, 27 1 show that 
condensed, superfluid clouds and normal clouds have dif- 
ferent frequencies of oscillation, leading to relative mo- 
tion between both components of the fluid. Finally, in 
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FIG. 1: (a) Experimental configuration to nucleate vortices in 
the BEG, showing the anti-Helmholtz coils slightly misaligned 
from the condensate long axis and (b) pictorial view of the 
expected movements/deformations suffered by the condensate 
cloud. 



Ref. [2J|, where we describe the vortex nucleation proto- 
col, we directly observe the quadrupolar modes as well 
as vortex structures in the interface of the dense core 
that represents the condensed cloud and the outer cloud 
(Fig.7 of Ref. [ill). We thus believe that this is the effec- 
tive mechanism contributing to nucleate the vortices and 
eventually lead to the turbulent state in the oscillating 
excitation. 

When the amplitude and/or excitation time is in- 
creased above some given values a completely different 
regime takes place. We observe a fast increase in the 
number of vortices followed by a proliferation of vor- 
tex lines not only in the original plane of the individual 
vortices, but distributed for the whole sample, covering 
many directions. An example of such a regime can be 
seen in Figl2^. In fact, equivalently to the already re- 
ported instability in the vortex formation pj], here we 
also observe an evolution from distinct vortices to the 
turbulent regime, where many vortices without preferred 
orientation are formed. Following the nomenclature of 
Tsubota et aL[§], we call this regime a vortices tangle 
configuration. Besides the fact that the image in figH^ 
comes from absorption in a single direction, careful obser- 
vation allows one to identify lines crossing the sample in 
many angles. For the purpose of visualization, together 
with the real absorption image we plot lines that indicate 
the presence and position of the vortices in the turbulent 
cloud in FiglJjD. 

The presence of quantized vortices non regularly dis- 
tributed along the whole sample in all directions is quite 
characteristic of turbulence in the quantum fluid and it 
is here taken as evidence for this regime. In fact, our 
way of nucleating vortices, where non-equivalent trans- 
lational and torsional movements of the quantum cloud 
occur in different planes of symmetry of the sample, is 
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FIG. 2: (a) Atomic optical density after 15ms of free expan- 
sion showing vortex structures spread all around the cloud 
resembling the vortex tangle regime proposed in Ref.^S] and 
(b) Schematic diagram showing the inferred distribution of 
vortices as obtained from image shown in (a). 

closely related to what is proposed in Q for combined 
rotations. 

An important question relies on which is the maximum 
number of vortices that can be introduced in a BEC in the 
turbulent regime. This number is certainly much smaller 
than in liquid Helium, mainly restricted by the available 
fluid volume (smaller in the BEC) and vortex core size 
(bigger in the BEC). Consequently some turbulence fea- 
tures might not be observable and/or be very different 
compared from that encountered in other cases of quan- 
tum turbulence, which involve large volumes containing 
many vortices, specially those related to the Komolgorov 
spectra, where large volumes are needed. Nevertheless, 
that low vortex density may allow one to make detailed 
observation of vortex reconnections and all the related 
phenomena such as cascade-like processes. 

Since the absorption imaging is destructive, in the 
present experiment it is not possible to observe the evo- 
lution of turbulence, vortex reconnections, Kelvin-wave 
cascades and all the dynamic behavior related to tur- 
bulence. Nevertheless, with the in-situ, non-destructive 
imaging technique of phase-contrast [32l |. it will be cer- 
tainly possible to observe the whole dynamics of the tur- 
bulent regime in cold atoms, from the transition from the 
regular vortex pattern to the turbulent pattern, going 
through vortex reconnections ending in the mechanisms 
that eventually lead to the decay of turbulence. Besides, 
imaging the atoms with this technique from two or more 
different directions, it is even possible to reconstruct [s^ 
the whole tridimensional vortex tangle distribution. 

Focusing on the turbulent regime, although we cannot 
predict the number of vortices formed nor their relative 
positions/orientations, we do not rely only on the visual 
signature of the cloud with plenty of vortices to assign 



the presence of turbulence. In fact, another remarkable 
feature is present in this cloud to confirm that it follows 
a completely different hydrodynamic regime: the sup- 
pression of the aspect ratio inversion during free expan- 
sion. This phenomena is only present after the turbulent 
regime takes place. The cloud does invert its aspect ratio 
when only a few regular vortices are present. 

As it is widely known, one of the most strong sig- 
natures of the quantum degeneracy in a gas of bosons 
is its asymmetric free expansion [3J]. In contrast to 
the isotropic expansion of a thermal cloud, an initially 
anisotropic trapped condensed cloud expands faster in 
the most confined directions than in the less confined, so 
a cigar-shaped cloud, when free to expand, turns rapidly 
into a pancake-shaped cloud with an inversion of its as- 
pect ratio. This phenomena is well explained theoreti- 
cally [1^ and is basically due to the domain of the inter- 
atomic interaction energy over the kinetic energy in the 
condensate. The former is proportional to the atomic 
density distribution which is anisotropic. 

In contrast to either the thermal and superfluid free ex- 
pansion behaviors, the turbulent cloud keeps its aspect 
ratio during the whole free expansion or, in other words, 
the aspect ratio inversion is suppressed. Still more re- 
markable is the fact that it does not evolve to the uni- 
tary aspect ratio as in a thermal gas, evidencing that this 
cloud is certainly not thermal, keeping somehow a quan- 
tum nature as it should be expected in a true quantum 
turbulent regime. 

Figure [3^ shows side by side several normal BEC 
clouds and turbulent clouds at different times of free ex- 
pansion, showing the described behavior. FiglSb shows 
the quantitative aspect ratio for the turbulent regime 
(black squares) and for the normal BEC (red circles). 
We believe that this remarkable characteristic is a new 
effect for atomic superfluids, possibly as a signature of 
the emergence of QT in this system. 

The change in the aspect ratio during free expansion of 
a vortex populated quantum fluid is a topic still not well 
explored in the literature. On the other hand, studies of 
hydrodynamic expansion of a rotating vortex free super- 
fluid cloud have already been presented [3^1 . In that case, 
it was observed that conservation of angular momentum 
combined to the irrotational flow, causes the expansion 
to take place in a distinctively way when compared to 
the non rotating superfluid. In our case, in contrast, the 
expanding cloud is populated of many vortices in many 
directions and just keeps its aspect. We believe that the 
presence of this random flow field along the condensate 
somehow acts preventing the known expansion regimes. 
A full theoretical understanding of this topic is still miss- 
ing and it would certainly shine some light on the true 
nature played by the flow of matter around the vortex 
cores and their interplay with interatomic interactions, 
with consequences on the hydrodynamics of the quantum 
fluid. Also, a full theoretical understand of this matter 
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FIG. 3: (a) Comparison between the free expansion of a nor- 
mal BEC showing the well known aspect ratio inversion (left) 
and a turbulent cloud (right) that keeps its aspect ratio while 
expanding, (b) Quantitative plot of the aspect ratios of the 
normal BEC cloud (circles) and the quantum turbulent cloud 
(squares). The lines are only a guide for the eye to stress the 
observed behavior. 



might help in the understanding of turbulence itself. A 
comparison of the expanding BEC and turbulent BEC 
(Figl3^) shows that the latter is typically much bigger 
for longer expanding times, which may be related to the 
decay of the turbulence through vortex line reconnection 
and excitation of Kelvin-waves and possible existence of 
cascades phenomena [37[ , leading to a heating of the sam- 
ple. 

In conclusion we have shown experimentally the emer- 
gence of quantum turbulence in an atomic quantum fluid 
evidenced by the presence of vortex cores as well as vor- 
tex lines in the absorption imaging non regularly dis- 
tributed across the sample. These configurations are 
known as vortex tangles, in contrast to the crystalliza- 
tion regimes, and are an important characteristic of the 
turbulent regime. Also, we have shown that the turbulent 
cloud does not expand like either a quantum or a ther- 
mal gas, but keeps a fixed aspect ratio. This expansion 
regime has never been observed and is also believed to 
be a signature of the presence of QT. With the assump- 
tion that the absence of inversion is a feature indicating 



turbulence in the quantum fluid it could also be used to 
indirect quantify this transition as well as being an in- 
dicative of a completely new hydrodynamic regime of the 
atomic superfluid. 

Experiments to observe QT in BEC are being carried 
out also at NIST, in the group of K. Helmerson [38| . 
We acknowledge helpful suggestions from J. Dalibard, 
M. Modugno, W. Philipps and C. Salomon and technical 
support of R.F. Shiozaki. 
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